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It is known that it is possible to encode a logical qubit over many physical qubits such that it 
is immune to the effects of collective decoherence, and it is possible to perform universal quantum 
computation using these 'decoherence-free' qubits. However, current proposed methods of perform- 
ing gates on these encoded qubits could be difficult to implement, or could take too much time to 
perform. Here we investigate whether exploiting ring-exchange interactions, which may be naturally 
present, can simplify the implementation of these gates in any way. Using a ring exchange inter- 
action, we have found a way to create a controlled-Z gate on the 4-site decoherence-free subspace 
and the 3-site decoherence-free subsystem, with high fidelity. The 4-site gate has interactions which 
are simpler than existing proposals, and the 3-site gate is also feasible although not necessarily an 
improvement on existing 3-site gates. We also investigate how timing errors and magnetic field 
fluctuations affect the fidelity of these gates. 
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I. INTRODUCTION 

Quantum computers offer a speed up compared to con- 
ventional classical computers, in problems such as factor- 
ing numbers flj or searching an instructured database [2] . 
One of the biggest obstacles to realising a useful, scalable 
quantum computer is making it robust against noise and 
errors, so that the coherent quantum information can be 
long-lived. There are two approaches to this problem: we 
could take a 'software' approach, where we encode our 
logical qubit over several 'physical' qubits in such a way 
that errors can be spotted and corrected easily; exam- 
ples of this include the nine-qubit Shor code [T] and the 
7-qubit Steane code L 3 S . The other approach is the 'hard- 
ware' approach, in which the logical qubit is encoded over 
several physical qubits in such a way that the noise has 
minimal effect in the first place. One particular example 
of this is where we encode the logical information in part 
of the full Hilbert space in which the noise has no effect; 
such a subspace is called a 'decoherence-free subspace' 
(DF subspace) or more generally a 'decoherence-free sub- 
system' (DF subsystem) [3HZ]- 

It was shown in [BHS] that universal quantum compu- 
tation can be performed inside a DF subspace or sub- 
system. To do this, we must be able to perform arbi- 
trary single qubit rotations as well as gates between two 
logical qubits (such as a controlled-Z gate) [§]• For en- 
vironmental noise which acts uniformly over all physical 
qubits ('collective decoherence'), explicit gate sequences 
which would realise this have been found in [5J [TOHT5] for 
qubits encoded over 3 or 4 physical qubits (forming a 3- 
site DF subsystem and 4-site DF subspace, respectively). 
Performing single qubit rotations in the 4-site DF sub- 
space and 3-site DF subsystem is relatively straightfor- 
ward (see e.g. [HHH]), however creating two-qubit gates 
is not; the methods found so far to perform two-qubit 
gates involve either interactions that could be difficult to 



create in an experiment [H [TT] , or involve large numbers 
of gates to be switched on and off sequentially (e.g. 50 
exchange gates in 27 time steps in [T^] for the 4-site DF 
subspace, and 22 gates in 13 time steps for the 3-site 
DF subsystem in |13j). or use perturbative/complicated 
control sequences to create these gates [14] . 

Our aim is to simplify the existing methods of univer- 
sal quantum computation in a 4-site DF subspace, by 
finding alternative ways to perform 2-qubit gates using 
less operations, and using operations which are easier to 
implement (i.e. require less control), so that realisation of 
an encoded qubit in a 3-site DF subsystem or 4-site DF 
subspace might be more attainable in an experiment. We 
are also interested in seeing whether or not the presence 
of ring exchange interactions can lead to simpler gates on 
these encoded qubits. This work is partly motivated by 
the recent experimental advances in realising quadruple 
quantum dots in a square configuration |15j . 

The paper is laid out as follows: In sec.llllwe introduce 
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some important background material. In sec. Ill we dis- 
cuss our main result constructing two qubit gates in the 
3-site DF subsystem and 4-site DF subspace, and then 
in sec. [V] we test the performance of these gates when 
timing errors or random magnetic fields are present. 



II. BACKGROUND INFORMATION 

A. Decoherence-free subspaces and subsystems 

Consider encoding a logical qubit over a number of 
physical qubits. If each qubit interacts identically with 
the environment (by being close enough together) we 
call the resulting decoherence collective decoherence. In 
this paper we will focus on DF subspaces and subsys- 
tems which are resistant to collective decoherence only. 
There are two types of collective decoherence we consider; 
the first is collective dephasing, for which the system- 
environment interactions do not change population of en- 
ergy levels, but change the coherence between them, and 
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for which the master equation takes the form [H)J E] : 
dp 



dt 



T([S zP ,S z ] + [S z ,pS z }), 



(1) 



where S z :— J2 n °« ^ s * ne collective z-operator, and 
ci?y,z are the x, y, z Pauli spin matrices acting on the 
n th physical qubit. The second type of decoherence we 
consider is collective amplitude damping, in which popu- 
lation transfer can occur. Then the master equation for 
this process takes the form [16] : 



dp 

dt 



it [S+S-, P ] + a ([Sip, S+] + [S-,pS+]) , (2) 
where S± = \{S X ± iS y ), and S x := 2„cri n) , 



S,, 



J2 n a y ■ The condition for states \j) to form a DF sub- 
space T>, given a master equation with Lindblad opera- 
tors L a is [6]: 



L a \j) = c a \j) for all \j) e V. 



(3) 



If the collective spin operator is S 2 = (S x ) 2 + (S y ) 2 + (S z ) 2 
with eigenvalues ms(ms + l)h 2 , and eigenvalues of S z 
are m z h, then to satisfy eqn. [3] for collective dephasing 
we need degenerate states which have the same ms value 
and the same value of m z |18j . To satisfy equation [3] for 
collective amplitude damping, we need states which an- 
nihilate when acted on with the SI and S+S_ operators. 
This is satisfied for states with m z = —ms |19) . 

The 3-site encoding which we use in this paper acts 
as a DF subsystem for collective dephasing but not for 
collective amplitude damping. It is defined over four or- 
thogonal states, defined by the eigenstates of the three 
operators S 2 , S z and S 2 2 , where S? • is the total spin 
operator acting only on physical qubits i and j, and has 
eigenvalues mij(mij + l)h 2 . Then in terms of the quan- 
tum numbers mg, m, mi^, the states making up the 
3-site DF subsystem are: 

\m 8 = 1/2, m, - 1/2,7^,2 = 0) = |<F)i 2 |0) 3 (4) 
11/2,-1/2,0) = k/r)i 2 |l) 3 (5) 

|l/2, 1/2, 1) = -^(72|T + ) 12 |1) 3 - |T ) 12 |1) 3 ) (6) 

|l/2, -1/2, 1) = -^(|r ) ia |l) 3 - V2|T_) 13 |0) 3 ). (7) 

Here we define the singlet states on qubits i and j as 
■— (|01)jj — |10)jj)/\/2 and the triplet states 



\T+)i 



|0Q)y, |T_) t3 



11 



To) 



(101), 



\10)ij)/y/2. We define the logical zero \Q)l in this 3-site 
subsystem to be an arbitrary combination of the first two 
states, i.e. |0) L := 5|l/2, 1/2,0) + -y|l/2, -1/2, 0), whilst 
the logical one is a superposition of the last two 
states \0} L := 5| 1/2, 1/2,1) + 7|l/2, -1/2, 1). The ac- 
tion of collective dephasing on this encoding can change 
the values of S and 7, but it will not couple states with 
different values of mi. 2, so the information is preserved. 



The 4-site encoding we use in this paper acts as a DF 
subspace for collective dephasing and also for collective 
amplitude damping. It is defined over the two states 
\ms = 0,m z — Q,mi2 = TO3.4 = 0) and \m s = 0,m z = 
0,mi 2 = m 3 A = 1): 



|0l) := 
■= 



12\lp )m 

\T + ) 12 \T\ 3i - |To) 12 |ro) 34 



1 

T-)i 2 \T+) 34 



(8) 



The 4-site encoding has one additional advantage over 
the 3-site encoding; it also functions as a supercoher- 
ent qubit [50]. Supercoherence would allow resistance to 
single qubit errors in any direction, and the mechanism 
for achieving this is as follows: When single qubit errors 
along any direction are applied to the ms = states in 
eqn. [HJ it is accompanied by a change in the ms value 
by 1 ^20) . In order to use this to create a supercoherent 
qubit, we switch on the Hamiltonian Hsc'- 



Use = Jsc ^ Ei : 



(9) 



where Eij := afaj + of <rj + of crj, and the sum is over 
all pairs of the 4 physical qubits. With this Hamiltonian 
switched on, the ms — states are degenerate and lowest 
in energy, with an energy gap between the 7715 = states 
any other states. Thus any decoherence process acting 
on individual physical qubits involves a change in energy 
of the encoded qubit, and will lead to a transfer of en- 
ergy from the environment to the system, which we can 
inhibit by cooling the environment. Thus supercoherent 
qubits would be very useful as quantum memories, and 
it was argued in |20j that computation with supercoher- 
ent qubits could be performed provided the interaction 
strength between qubits was small enough compared to 
Jsc (leading to a trade-off between the speed of opera- 
tions and the robustness against errors). In this paper, 
we will not aim to make our interactions supercoherent 
as well (i.e. we envisage a protocol in which we use the 
supercoherent mechanism as a means to reliably store 
information, but turn off the supercoherent Hamiltonian 
Hsc when we interact endcoded qubits). 

For this investigation, we consider physical qubits ar- 
ranged in a regular formation: three spins in an equilat- 
eral triangle for the 3-site DF subsystem, and four spins 
in a square plaquette for the 4-site DF subspace. The in- 
teractions we consider for both encodings are interactions 
between the middle 4 spins (see fig. [T for an illustration 
of this) . As an example of why we might want to simplify 
the existing two-qubit gates in the 4-site DF subspace, 
see fig. [2] for an illustration of one such interaction used 
in [SI [TT] to create a gate. The Hamiltonian for this in- 
teraction is: 



H = 3E- 



34 



E 



13 



-^32)- 



(10) 
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Given the large difference in couplings between qubits 3 
and 4 compared to the other qubits, and the absence of 
coupling between qubits 1 and 4, this could be a chal- 
lenging gate to realise (note that although this only acts 
within one encoded qubit, this Hamiltonian is turned on 
whilst we are out of the logical subspace, and so is not 
simply a local unitary operation). 



4-site encoding 




Logical Qubit 1 Logical Qubit 2 



3-site encoding 




Logical qubit 1 Logical qubit 2 



FIG. 1. An illustration of the two geometries of qubtis we 
consider when constructing our logical qubits. Filled circles 
represent physical qubits, and solid lines illustrate the kind of 
interactions we consider when looking for a gate in this paper. 
The top diagram shows the layout for the 4-site encoding and 
the bottom diagram shows the layout for the 3-site encoding. 




Qubit 1 



FIG. 2. An illustration of interactions used in [11] to inter- 
act two 4-site encoded qubits, which could be challenging to 
implement a 4-site encoding in an experiment. Filled cir- 
cles represent physical qubits, and solid lines represent one 
of the interactions used in. Note that this occurs when the 
qubits are taken out of the logical subspace into the larger 
14-dimensional singlet subspace over 8 qubits, and so is not 
simply a local unitary transformation. 



B. Ring Exchange Interactions 



When constructing our gate, we will include 'ring 
exchange' interactions. These interactions, which are 
used to explain e.g. excitations in La-iCuO^ and 
become important in electrons forming a wigner crys- 
tal [23l |24] , appear as corrections in the exchange Hamil- 
tonian due to higher order hopping processes between 
different sites. They have also been investigated in the 
context of quantum computing [25'.427j . and it is clear 
from these papers that ring exchange processes should 
not be ignored. Including ring exchange terms, the mod- 
ified Hamiltonian for 4 spin- 1/2 particles located at sites 
1,2,3 and 4 becomes: 

3 2 
H = JnN ^ E ntn+ i + JnNN ^ £n,n+2 
n—1 n—1 

+ Cl234[-El2-E'34 + E14E23 — E13E24] 
+ Cl324[-El3-E'24 + E14E23 — E12E34] 
+ Cl342[-El3£'24 + E12E34 — E14E23] 

+ (ID 



where t is the tunneling coefficient, U is the on-site 
Couloumb interaction, Jij is the exchange coupling be- 
tween sites i and j, and Eij is the exchange interac- 
tion defined in sec. |H] The ring exchange coefficients are 
Cijki oc tijtjktkitu where tij is the tunneling coefficient 
between electron sites i and j (which for simplicity takes 
values t or 0), and non of the indices i,j,k,l are equal. 
This means these ring exchange terms will only appear 
when there are exchange terms present which form a loop 
(i.e. when a physical qubit is indirectly coupled to itself). 
For example, if we have four electrons with exchange in- 
teractions between electrons 1-2, 2-3, 3-4, and 4-1, then 
there will be ring exchange terms present. In [25] . they 
also find that the presence of magnetic fields changes 
the coefficients Cijki and introduces three-body terms, 
with couplings that depend on the magnetic flux passing 
through 3 or 4-site loops. Here we assume that these 
magnetic fields are low enough so that the magnetic flux 
has negligible impact and three body terms can be ig- 
nored, and we leave the effects of larger magnetic fields 
to later work. 

With four qubits arranged on a plaquette, it possi- 
ble to make all of the couplings uniform (i.e. Jnn — 
J nnn = J), and since this will simplify things consider- 
ably, whenever we need to include ring-exchange terms 
we will assume that this is the case. Taking uniform 
interactions in eqn. [IT] results in a 'symmetric' version 
of the Hamiltonian, Hs, which takes the form derived 
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H S = J(H D +H X ) + J H 

4 2 
= J J]] E n n+ i + J E n n+2 
n—1 n—1 

+ Jo[E 12 E 3 4 + E14E23 + B13-B24] 
= J(H n +H x +aH ), 



(12) 



where H\j and ff x represent the nearest-neighbour and 
next-nearest-neighbour interactions, respectively, and 
Hq represents the ring exchange terms. In the final line 
we have defined a := Jq/Jnn = Jo/ J as the ratio of 
ring exchange terms to the nearest/next-nearest neigh- 
bour terms. Restricting ourselves to symmetric Hamil- 
tonians of this form simplifies things considerably; 
commutes with if □ and ff x , and as argued in |26] , this 
form of Hamiltonian contains enough degrees of freedom 
to fix all of the eigenvalues, and so we do not need to take 
into account any higher order terms, unlike the pertur- 
bative expansion in eq. 
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Since single qubit rotations of 
the encoded qubit can be performed by exchange inter- 
actions between physical qubits, creating a two-qubit in- 
teraction comes down to simulating a 4-body interaction 
of the form EijEki, and since the ring exchange interac- 
tions contain terms similar in nature to these, it seems 
plausible that we can use these to simplify the 2-qubit 
gates acting on the encoded qubits. 

We also introduce notation for two other Hamiltonians, 
which only have couplings in the vertical or horizontal 
directions: 



fix 
if II 

so that H a =H~+ Hu 



E\2 

Eu 



E34 

^23) 



(13) 
(14) 



III. TWO-QUBIT GATES 

We now look at creating a controlled- iT gate (CZ) be- 
tween two encoded qubits, which, along with arbitrary 
single-qubit rotations, will enable us to perform universal 
quantum computation [5]. Performing two qubit gates is, 
predictably, not as straightforward as single gates, mainly 
because we move into a much larger Hilbert space as 
soon as interactions between the two qubits are turned 
on. Since [S 2 ,Eij] — (provided all i,j are both qubits 
that S 2 operates on), if we start with both qubits in a sin- 
glet subspace, then when we interact them together using 
exchange interactions they will still form singlets over- 
all, but they will form singlets over two encoded qubits 
rather than just one, so we can limit ourselves to this 
singlet subspace to speed up our calculations. 

We also constrain ourselves to couplings which are re- 
alistically possible: we only couple sites which can real- 
istically be placed near each other (by e.g. placing two 
encoded qubits side-by-side, see fig. [I]) . The method we 



used to search for a quantum gate, once we had chosen a 
certain set of interactions to use, uses the invariant quan- 
tities found by Makhlin |28j . In this paper, two invariant 
quantities mi and are derived for 2-qubit operations. 
Given a 4 x 4 matrix M, we first transform M into the 
Bell basis, M -» M B = Q^MQ where 



Q = 



V2 



( 1 i 
i 1 
i -1 

V 1 -i 



(15) 



Then defining m — M^Mb, rai and 7712 are given by: 



mi(M) = (tr m) 2 detM t 
m 2 (M) = tr(m 2 ) det . 



(16) 
(17) 



We can find these two invariant values for any operator, 
and if they both match with the value for a CZ gate 
then the two gates are equivalent apart from some local 
operations. 



: 11 ;i : 



Logical Qubit 1 Logical Qubit 2 Logical Qubit 1 Logical Qubit 2 



Hn 



Logical Qubit 1 Logical Qubit 2 Logical Qubit 1 Logical Qubit 2 

FIG. 3. An illustration of the interactions we use to construct 
the gates, excluding ring exchange interactions which are diffi- 
cult to represent in this form. Filled circles represent physical 
qubits, and solid lines illustrate exchange interactions. 

The set of Hamiltonians that we used to construct a 
gate with were: 



{ fix , if || , fix , if □ , ff(3 } 



(18) 



where each of these Hamiltonians is defined in sec. Ill Bl 
and illustrated in fig. [3] (except for Hq which is hard to 
illustrate in this way). For each Hamiltonian in this set 
we define a corresponding unitary operator: 



exp 



. i n fin 7~ n 



exp (—iH n 6 n 



(19) 



e.g. Cf x = exp(— zfi x (9 x ), etc. 



where 6 n := J„T n /h. Note that all of the 8 n are phases, 
so we are free to add multiples of 2tt without changing 
the properties of the gate. We then multiply all of these 
unitaries together in some order to get the full gate op- 
erator: 



Y[U n = Y[exp(-iH n O n ). 



(20) 
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Before comparing this to a CZ gate, we first project into 
the logical subspaces of the 3- or 4-site encodings, so 
Utot —> Ulot — PjJJtotPh-, where Pl is a projector into the 
logical subspace (note that, in the 3-site DF subsystem 
case, subtleties arise as we must project from two states 
into one, i.e. states 1 1/2, ±1/2, 0) are projected into 
|0)l and states 1 1/2, ±1/2, 1) are projected into |1)l)- 
We then find the Makhlin invariants mi(U' tot ),m2{U' tot ) 
of the resulting unitary U' tot , and compare these to the 
Makhlin invariants m%(CZ), m^iCZ) of an ideal CZ gate 
using the following function 



/ m = ^K(CZ)-m i (U t ' ot )| 



(21) 



if we have within our gate sequential applications of the 
Hamiltonians H x , Hn and Hq with parameters 9 X , 
9q respectively, we can express this as the Hamiltonian in 
cqn. [12] plus an additional next-nearest-neighbour term 
U' x with parameter 9 X out in front (if 0' x is negative we 
can just add 2-7T since it is just a phase). Also note that, 
since we have set tjj = t@ , then a = 6q/9jj. This identity 
means that for the purposes of the minimisation we can 
keep the H x , Ha and Hq as if they were separate inter- 
actions with independent parameters and then combine 
them together at the end using the identity in eqn. |23| 
and when we do combine them the value of a is set by 
Bq/Ou. 



which gives us a measure of how close we are to a CZ 
gate, excluding local unitary rotations. Minimising over 
f m will give us possible ways to implement a gate (or will 
tell us if it isn't possible for the type of interactions we are 
considering). The local operations required to transform 
our result to a CZ gate will be easy to find compared to 
the difficulty of minimising over f rn , and so in this paper 
we focus on finding sequences which minimise f m . We 
also need to consider how far out of the logical subspace 
we are, so we define the leakage parameter as: 



L:=<i-Y \(i\U'\j}\ 



(22) 



Where \i), \j) € C denotes states in the four-dimensional 
logical subspace of the two encoded qubits, C — 
{|0l0l), |0 l 1l), |li0i), |1l1l)}- Using these measures, 
we implemented a genetic algorithm (see e.g. [29]), fol- 
lowed by a Nelder-Mead simplex search [3D] once we had 
narrowed our search down to a sufficient level. The same 
method has been used in [10] and [T2j . 

Finally we note an identity which makes this search 
over parameters easier to make. As noted in section [Tl B| 
we consider Hq such that it commutes with H x and H\j , 
which allows us to rearrange U x U\jUq as: 



U x U a U = e 



: e -iH n 9 n e -iH e Q 



= e -i(H x e x +H n 8 n +H e ) 

= e ~^ H x- ( x -0n)) e (-*ffx e a +H a e a +H e ) 

= e ~*(H x (0 X -8 D )) e -iJr D (H x +H n + (J /J)H ) 
= e -i{H x (0 x -e a )) e -i6a{H x +H a + a H o ) 



where 



= u' x u s . 



U = e -i[H n +H x +aH ] = e -iH s 8 s 



:-»□)) = e -i(. H x 



(23) 



(24) 



In the above we made tq = tq so that H\j and Hq op- 
erate at the same time, and recall that a — Jq/J and 
J\j = J x — J in order to use the symmetric form of the 
ring exchange interaction, Hq. Then using this identity, 



A. Constructing the gate 

We tried many combinations of the Hamiltonians given 
above, performing a genetic search for each one and then 
a Nelder-Mead search if the genetic search yielded f m ~ 
0.1. The following combination was successful for both 
the 3-site and 4-site encodings: 



U gate = U&U^UxUuUoU^Ug 

= u&u<PuLu a uV>i& 



(25) 



with corresponding parameter set: 

{0} - {0g\Bp,6 x , 9 U , 9 ,9f\ 9^}, (26) 

where the superscripts (1) and (2) are used to differen- 
tiate between the interactions used at the beginning and 
at the end of the gate. 

Using this combination of interactions, with parame- 
ters {9}3 for the 3-site qubit and {9}4 for the 4-site qubit 
(given in table |l|, we were able to make gates with a pro- 
cess fidelity of 1 — 1 x 10~ 5 and 1 — 3 x 10~ 6 respectively 
(for definition of process fidelity see sec.|V]). If we assume 
that for each interaction the coupling strength is limited 
to some maximum possible value J ma x , then we can find 
the total gate time T in units of h/J maxi as an indica- 
tor of how long this gate would take compared to other 
gates. Note that we do not simply add the parameters 
in table [I] since we apply the identity in eqn. [23] first. 
This gives gate times (in arbitrary units) of 17.7 for {9}^ 
and 18.4 for {9}^ (shown in table|l]in the row labelled T). 
This means we would expect our 4-site DF subspace gate 
to be faster than the gate in [12] (which has a total gate 
time of 52.1) but about as fast as the gate in [TT] (which 
has a total gate time of 18.7). On the other hand our 3- 
site DF subsystem gate takes longer than the gate in [13] 
(which has a total time of 9.9 in these units), especially 
if we make the ring exchange coupling more realistic (see 
sec. Evil. 
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TABLE I. Possible Combinations 





{Oh 


{0U 




5.4913836 


4.8245730 


ep 


0.4058571 


1.9550064 


e x 


1.5033263 


2.3010055 


On 


3.4556824 


5.0837756 


9q 


0.8481379 


0.8471519 


q( 2 ) 

II 






4 2) 


3.5520005 


1.8441184 


fm 


5 x 10~ 4 


3 x 10" 5 


L 


5 x 10~ 5 


3 x 10~ 5 


l-F p 


1 x 10 -5 


3 x 10~ 6 


T 


17.65 


18.37 


a 


0.245 


0.167 


a« 


0.087 


0.075 



IV. REMOVING SOME CONSTRAINTS 



The results in table U are limited in that the rela- 
tive sizes of the nearest-neighbour couplings, J and ring- 
exchange coupling J(j , are set by the ratio a. = JqJ J = 
9q /9a , since as we have seen in eqn. 23 we end up turn- 
ing on the Hamiltonian H x + Ha + aH^ for some time 
t\3 = 0\jh/ J- This means we are constrained to situations 
where we can set a to 9^ /6a, and applying this to {9}3 
and {9}i gives — 0.245 and a± = 0.167 respectively. 
We would expect a < 0.17 (see e.g. [ZOEiSO!]), so only 
q;4 seems like a reasonable value, and is too big to be 
realistically obtainable. However, since all of the 9 values 
are just phases, we are free to add factors of 2tt to any 
of them, so we can add multiples of 2ir to 6a in order 
to decrease the value of a in our implementation. Since 
9a = Jta/h, the trade-off is that increasing 6a corre- 
sponds to either increasing J or increasing ta , and under 
the assumption that we are using the strongest couplings 
possible, this really means increasing the gate time by 
2ith/J every time we add a multiple of 2ir. This is not 
ideal, but at least gives us a way of making a 3-qubit gate 
with a realistic value of a. Values of a in which multiples 
of 2ir have been added to &a are given in table |TJ and we 
label them as = 9 /9^\ where 6$ :=(9a + 2rwr). 
In order to make {9}3 more realistic, we have to add 2tt 
to 9a, giving o/ 1 ) = 0.087, which would mean a slight in- 
crease in gate time from 17.7 to 23.9 (a relative increase 
of 35%). Since {9}^ already has a low value of a, the 4- 
site gate is slightly more useful, as it can be implemented 
in less time. 



Clearly this has given us more flexibility, but wc still 
seem to be tied down to a few precise values of a^ and 
a±. To get around this, we notice that we can split up 



the ring exchange into two parts: 
U s = exp (-iH x 9 a ) exp {-iH a 9a) exp (-iH 9 ) 
= exp (-iH x {9 D + 9 h a )) exp (-tIT n (0& + 9 b D )) 
exp(-iH {9 a Q +9 b )) 

= exp (-iH x 9 a ) exp (-iH a 9 a ) exp (-iH 9^) 
exp (-iH x 9 a ) exp (-iH a 9 a ) exp (-iH 9 b ) 
= U a x U^U b x U D U b , 
■= U a s U b s 



(27) 



0q . 



(28) 
(29) 



so we now have two ring exchange interactions, Ug and 
Ug which have the same form as H$ but with different 
values of a. Now since 9 n := J n t n /K, and since UgU s — 
Us and all of the terms commute, this means that: 

J a t a + J b t(, 

where 9g, 6q are the parameters we found in the search 
in sec. Illll and we have defined 9 a = J a t a , 9 a = J b tt,. 
For simplicity, we take J a = J b = J, without loss of 
generality, since we are free to scale these parameters as 
we wish, provided we scale the corresponding Jq values 
correctly. Since 6a and 9q are phases, we are free to add 
multiples of 2?t to these values, and so using the definition 

of 6^ ] fr om above, we can write: 



n(n) 



J a t a + J b tb 



Rearranging these equations gives us: 

at") 



tn 



h = 



a(n) 
□ 

J 

n(n) 

J 



a a — a b 

a a - a (n) 
a a - ab 



(30) 
(31) 

(32) 
(33) 



where a a '■= Jq/J, «6 := Jq/J, an d is defined 
above. For t a and tb to be positive, we need couplings 
such that a a > > ab- So this tells us that if wc 

are able to control the relative strengths of the ring and 
nearest-neighbour terms, and if we could get them such 
that a a > oS- n ' > ctb is satisfied, then regardless of what 
the actual values of a a and ab are, we can create the 
C Z gate (at the expense of adding more interactions and 
thus increasing the time of the gate) . 

For the remainder of this paper we will assume that we 
have made the replacement 6a — > 6a + 2-7T for the 3-site 
qubit, so the corresponding value of a is more realistic. 
We make no change to the value of 9a for the 4-site qubit, 
since it already has a realistic value of a. 



V. PERFORMANCE 

We now look at the performance of this gate under the 
influence of noise. The two types of noise we consider are 
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errors in coupling strengths or timing of the gates (i.e. 
random errors in the {9} parameters when implement- 
ing each of the Hamiltonians in cqn. 25 ) and fluctuations 
in magnetic fields acting on the qubits during the gate 
implementation (which would normally be protected by 
the decoherent-free or supercoherent properties, but are 
not whilst our gates are being implemented). The rea- 
sons for picking these particular types of errors are that 
one of the most promising systems in which to imple- 
ment these encoded qubits is in arrays of quantum dots, 
and there have been several advances towards achieving 
these, e.g. [TU [3T]. These quantum dots are suscepti- 
ble to errors in exchange coupling due to charge fluctua- 
tion [32] (which we are modelling by fluctuations in the 
{9} values) and fluctuations in external magnetic field 
due to the nuclear spin bath or stray magnetic fields (see 
e.g. [33]). When looking at these magnetic fluctuation 
errors, we look at two cases; in the first case we assume 
that magnetic field fluctuations are roughly uniform over 
each encoded qubit (so that there is collective decoher- 
ence acting on it, see sec.|l), but the magnetic fields act- 
ing on different encoded qubits are different in magnitude 
and direction. In the second case, we consider magnetic 
field fluctuations acting independently on each physical 
qubit (i.e. a situation where a supercoherent qubit would 
be more appropriate) . For all errors, we assume that the 
time scale for the fluctuations is large compared to the 
time to perform the gate, which is typically the case [33J- 

To measure the size of the errors, we calculate the gate 
distance using the techniques of quantum process tomog- 
raphy [3(3 [37] ■ Suppose we have a process acting on a 
state p in a d-dimensional Hilbert space, such that p is 
mapped to £ (p). We can describe this process as: 



£{ P ) = Y J EwE\ 



(34) 



where {Ei} are Kraus operators which satisfy J^i EiE\ < 
1, with equality iff £ is a trace-preserving map. For our 
process, we consider performing our sequence of oper- 
ations which are locally equivalent to a CZ gate, then 
performing some local operations on each qubit, and 
then performing a measurement on each of the encoded 
qubits to determine whether the encoded qubit is in the 
logical basis or not. The practicalities of performing 
this measurement is not covered in this analysis, since 
the main reason we use it is to simulate how our gate 
changes with noise (and projecting our gate into the 4- 
dimensional logical subspace allows us to easily compare 
our gate to the 4-dimensional CZ gate). Similarly we 
expect the local operations to be easy to find [TU], so 
we do not find these explicitly; instead we take the gate 
with zero noise U to t corresponding to one of the parame- 
ter sets {9}, and find the matrix V which diagonalises it 
(i.e. e^V^UtotV — CZ, where $ accounts for any phase 
terms). Then once we have a gate with noise added, U' , 
we still use V and $ to convert U' into a noisy CZ gate 
(i.e. e^V^U'V « CZ), since we are restricted to always 



using the same local operations (i.e. we cannot change 
our local operations since we don't know what the noise is 
doing to our system). These local operations V and pro- 
jection onto the logical subspace combine to give a single 
Kraus operator E = -^Pie^V^U'V for each dataset 
{9}, where Pl is a projection into the two-qubit logical 
subspace and Af normalises the density matrix. To find 
the process fidelity, we decompose this Kraus operator as 
E = ~}^ m a m A m , where {A m } is a set of 4 x 4 matrices 
that form an orthogonal basis under the Hilbert-Schmidt 
inner product, i.e. tr(A\ n A n ) = 4S mn . Then we can 
rewrite eqn. [34] as: 

£(p) = EpE^ 



— ~Jj J]] XmnAmpA^ 



(35) 



where Xr. 



a m a* n is the 'process matrix'. In order to 



find the normalising factor, we note that tr(£(p)) = ^p-, 
so that J\f — tr(x)- Then we can absorb this into % by 
letting xAKx) ~^ X- Once have our operation in terms 
of this normalised process matrix x, we can then compare 
it to the process matrix of an ideal CZ gate, xcz, using 
the process fidelity F p [37]: 



F p (x,Xcz) = tr (^J ^/xczxVxcz^j 



(36) 



The process fidelity is a natural choice for measuring the 
accuracy of a gate, since 1 — F p has the interpretation as 
the upper bound of the average failure probability p e [36j . 



As noted in sec. IV the gate sequence for the 3-site qubit 
only gives a realistic value of a if we add an extra 27r to 
and so for all of the following simulations we will 
make the transformation — > + 2n for the 3-site 
qubit. 



A. Timing Errors 

For the first type of noise, we added random Gaussian 
noise to each of the parameters: 

6 n ->0 n + 8O n , (37) 

where 89 n is sampled from a gaussian distribution with 
mean and standard deviation e. This simulates ran- 
dom errors in values of the couplings between sites or 
timing when implementing gates. An example of the for- 
mer type of error is charge fluctuations in quantum dot 
systems (see e.g. [32]). Over the range e € [0,0.05] of 
e, we calculated the process fidelity over 150 iterations 
taken from a normal distribution with standard devia- 
tion e and mean 0, finding the average over all of these 
iterations (note that in this case, since all exchange op- 
erators commute with S 2 for any number of sites, we will 
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not leave the 14-dimensional singlet subspace in these 
calculations, which speeds up the process considerably). 
The results are shown in fig. |4j with only one set of re- 
sults shown since both the 3- and 4-site qubits give very 
similar results. The process fidelity falls off slowly and 
with stays above 0.99 for e < 0.02; a reasonable estimate 
of these fluctuations in gate couplings would be around 
0.01 |32j . at which point both gates have fidelity ~ 0.995, 
so we can see that these gates still have high fidelity even 
with this level of noise. 




FIG. 4. The average process fidelity F p when performing 
the gate with random fluctuations in gate times, where the 
fluctuations have standard deviation e. Only results from one 
data set, {0} 4, is shown, since both data sets showed a similar 
behaviour. A fourth-order polynomial fit to the data is also 
shown. 



B. Nuclear spin bath fluctuations 

If we were to implement this supercoherent qubit in 
a quantum dot, then there could be random magnetic 
field fluctuations due to the nuclear spins in the material 
or stray magnetic fields. We studied two scenarios which 
could occur; in one case we looked at the effects of having 
random field fluctuations which are uniform over all phys- 
ical qubits, but which may vary between encoded qubits 
(so that the collective decoherence assumption is valid for 
single encoded qubits, but when we interact two of these 
the assumption is not valid) . In the other case we consid- 
ered having independent random magnetic field fluctua- 
tions on each physical qubit. For both of these cases, we 
followed the arguments in [33J [35] , using the quasistatic 
approximation in which the magnetic field from the nu- 
clei B stays constant over the time we perform the gate, 
has random direction and has magnitude \B\ following a 
Gaussian probability distribution; 

P(\B\) = {27rB l /2 exp(-|B| 2 /2BL c ) (38) 

where B nuc is the standard deviation in fluctuations of 
magnetic field. We took the average over 150 iterations, 
with each iteration having a different magnitude and di- 
rection of magnetic field sampled fom the above Gaussian 



distribution. The magnetic field strength was taken rel- 
ative to the nearest-neighbour coupling strength J. In 
each case we found how the process fidelity varied for 
the 3-site and 4-site encoding. The results are shown 
in fig. [5] and fig. [6j Based on the exchange values given 
in [38] and the values for B nuc given in [33] . we would 
expect B nuc /J to be at most ~ 0.01. With this kind 
of nuclear field present, the 4-site encoding achieves a 
process fidelity of ~ 0.99 if there are errors across the 
encoded qubits, or ~ 0.98 if there are errors on each in- 
dividual physical qubit. The 3-site encoding achieves a 
process fidelity of ~ 0.97 if there are errors across en- 
coded qubits, or ~ 0.91 if there are errors on physical 
qubits. So overall, as we might expect, the qubits are 
much more robust to errors which are uniform on the 
physical qubits, but less robust to errors which vary be- 
tween physical qubits. Also, the gate on the 3-site qubit 
performs worse than the 4-site one, most likely because 
of the longer time required to perform the gate. Higher 
fidelities could be achieved by applying techniques to re- 
duce the effects of fluctuations in nuclear spin, such as 
the methods presented in [35] and [20] • Also a more ex- 
haustive search might yield parameters for the 3-site gate 
with lower gate times. 




FIG. 5. Average process fidelity F p as a function of B nuc / J, 
the magnetic field strength, with the fields pointing in random 
directions on each plaquette. The range is over B nuc /J = — > 
0.02. {6*}3 and {9}4 indicate results corresponding to a 3-site 
or 4-site qubit, respectively. A fourth-order polynomial fit to 
the data is also shown. 



VI. CONCLUSIONS 

We have demonstrated a simple way to implement a 
CZ gate in the 4-site DF subspace and the 3-site DF 
subsystem, using a sequence of 6 operations, excluding 
local operations, and including ring exchange interac- 
tions. The advantage of the 4-site gate is that we use 
more uniform Hamiltonians than in [11) . and the gate 
time is substantially shorter than in [12], which is es- 
pecially useful given the current advances in creating 4- 
electron quantum dots [T5]. This makes a 4-qubit DF 
subspace qubit much more attainable, and closes the gap 
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6 0.605 O.'Ol 0.015 0.'02 



FIG. 6. Averagee process fidelity F p as a function of B nuc /J, 
the magnetic field strength relative to the exchange coupling, 
with the fields pointing in random directions on each physical 
qubit. The range is over B„ uc /J = — ¥ 0.02. {6)3 and 
{9}4 indicate results corresponding to a 3-site or 4-site qubit, 
respectively. A fourth-order polynomial fit to the data is also 
shown. 



in practicality between the 4-qubit subspace and 3-qubit 
subsystem encodings. The 3-site gate we have found here 
does not appear to improve much on the exchange-only 
gate in [T3] , but further searches could yield a gate which 
takes less time. 

The process fidelity achieved was F p = 1 — 1 x 10~ 5 for 
the 3-site qubit and F p — 1 — 3x 10 -6 for the 4-site qubit, 
giving an upper bound on the average failure rate p e of 
1 x 10 -5 and 3 x 10 -6 respectively. We introduced errors 
when performing these gates, to simulate errors in cou- 
pling strength or gate times, and to simulate fluctuations 
in magnetic field due to some external environment, e.g. 



nuclear spins in a quantum dot. With nearest-neighbour 
exchange coupling J, we found that the 4-site gate kept 
F p above 0.99 even with nuclear fluctuations of around 
0.1% of J, or with timing errors of up to around 1% of 
J/h. This gate could be useful for current experiments 
in which there is limited control over the interactions in 
the qubit, or for which we would like to implement a 
CZ gate quickly, and with a smaller number of gates. In 
future work we would like to perform more exhaustive 
searches to find better fidelity gates with lower values of 
a, using the more general form of the ring exchange in- 
teraction rather than the one we have used here, and also 
including the effects of magnetic flux on the couplings as 
reported in (25]. From these results, it seems that includ- 
ing ring-exchange interactions may allow a simplification 
in the gate, which we might intuitively expect since the 
most direct CZ gate on these encoded qubits would in- 
volve a four-body interaction, which is present in the ring 
exchange terms. In future we would to investigate this 
more rigourously, to see if gates on these encoded qubits 
involving ring exchange terms outperform gates without 
ring exchange interactions. It would also be interesting 
to see if we could extend the techniques such as dynami- 
cal decoupling and leakage reduction, previously applied 
to 3-site encoded qubits in [13l |4T] , to this 4-site qubit 
gate. 
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